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The structural, electronic, and vibrational thermodynamic properties of the (ZnO)n (n=2-16) 
clusters are studied using density functional - full potential computations. The results show, small 
clusters up to n = 9 stabilize in the 2D ring shape geometries while the larger clusters prefer the 3D 
cage Uke structures. The ring to cage structural cross over in ZnO clusters is studied by investigating 
the behavior of the Zn-O-Zn bond angle, the Zn-0 bond strength, and the number of bonds in the 
systems. It is argued that 12 is the lowest magic number of ZnO clusters at ground state, while finite 
temperature vibrational excitations enhance the relative stability of the (ZnO)9 cluster and make 
it a magic system at temperatures above about 170 K. The obtained electronic structure of ZnO 
clusters before and after applying the many-body GW corrections evidence a size induced red shift 
originated from the ring to cage structural cross over in these systems. The behavior of the extremal 
points of electron density of the clusters along with the extrapolated cluster binding energies at very 
large sizes may be evidences for existence of a metastable structure for large ZnO nanostructures, 
different with the bulk ZnO structure. 
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I. INTRODUCTION 

In the family of the diluted magnetic semiconductors, 
the ZnO based alloys exhibit paramount technological 
advantages, because of their direct and wide band gaps, 
large exciton binding energies, large piezoelectric con- 
stants, strong luminescences, large non-linear optical co- 
efficients, high thermal conductivities, and availability of 
large single crystals [l|. The direct band gap of ZnO 
(3.37 eV i) is simificantly larger than the direct gap 
of GaAs (1.43 eV Q), while its exciton binding energy 
(60 meV) is much greater than the wide gap semiconduc- 
tor GaN exciton binding energy (25 meV)Q. 

Recently, ZnO nanostructures have attracted great at- 
tention for novel technological applications [1, Q • ZnO 
clusters have been synthesized by using the laser abla- 
tion and the electroporation of unilamellar vesicles [3] 
techniques and characterized by using the reflection time 
of flight (TOF) mass spectrometry. Pioneering theoret- 
ical studies on ZnO clusters were performed by Matx- 
ain et al. in the framework of density functional theory 
(DFT) with the hybrid B3LYP functionals 8] and time 
dependent DFT [|]. They found that the 2-7 units clus- 
ters stabilize in two dimensional ring shape geometries 
while the larger clusters favor three dimensional config- 
urations. Moreover, the 3D ZnO clusters were found 
to have smaller excitation energies compared to the 2D 
ring like clusters. Wang et aZ.[l^ studied ZnO clusters 
by using DFT with the generalized gradient functionals 
and concluded that the structural cross over from the 
2D ring geometry to the 3D cage/tube configurations 
occurs in the 8 units cluster. The (ZnO)i2 cluster has 
been reported as the most stable system among the 2-18 
units ZnO clusters. More recently, they proposed sodalite 
structure being made of coalesced (ZnO) 12 cages as a 
metastable geometry for larger ZnO nanostructures [ll[ . 

In this paper we apply DFT computations to investi- 
gate the structural, electronic and vibrational properties 



of the (ZnO)„(n=2-16) clusters. The employed computa- 
tional method are discussed in the next section. Then the 
structural behavior of ZnO clusters are studied in terms 
of atomic configurations, bond lengths, bond angles, ex- 
tremal points of electron density, and binding energies of 
the systems. Next two sections are devoted to the elec- 
tronic and vibrational properties of the most stable ZnO 
clusters. The obtained vibrational spectra is used to cal- 
culate some thermodynamic properties of the systems. 
Our conclusions are presented in the last section. 



II. METHOD 

The electronic structure calculations and geome- 
try optimizations of the present work were performed 
in the framework of density functional theory, us- 
ing the all-electron full-potential "ab initio molecular 
simulation" (FHI-aims) package [l^]- This code employs 
the numeric atom-centered orbital basis functions to 
achieve both computational efficiency and accuracy for 
investigation of the cluster as well as the periodic ge- 
ometries. The calculations were performed in the scalar 
relativistic limit by ignoring the spin-orbit interaction 
which is expected to be weak in ZnO materials [13]. We 
applied the BLYP-GGA functional composed of Becke 
exchange 14 1 and Lee- Yang-Parr correlation fl^ for 
exchange-correlation energy and the Broyden-Fletcher- 
Goldfarb-Shanno (BFGS) algorithm with a convergence 
criterion of lO^'^ au and force accuracy of 10""* Ry/bohr 
for cluster geometry relaxations. For calculating the 
vibrational frequencies and IR spectra the convergence 



criterion and force accuracy were increased to 10^'^ au 
and 10^^ Ry/bohr, respectively. The fundamental vibra- 
tional modes of the clusters were determined by calculat- 
ing and diagonalizing their dynamical matrix. In order to 
obtain reliable excitation spectra, the many body pertur- 
bation based GW correction was applied to the converged 



electronic structure within the BLYP-GGA functional. 

For more accurate description of electron density in 
ZnO clusters, we define topological critical points which 
are the extremal points of the scalar field of electron den- 
sity. By using the eigenvalues of the matrix of the second 
derivative of electron density (Hessian matrix) , four kinds 
of the critical points (CPs) are identified: the Nuclear 
CPs which are the local maxima that usually occur on 
nuclei, the Bond CPs which are the saddle points of elec- 
tron density between two neighboring atoms, the Ring 
CPs which are the second type of the saddle points ap- 
pearing in a ring configuration of atoms, and the Cage 
CPs which are the local minima of electron density occur- 
ring inside a cage configuration of atoms. The number of 
the critical points (/) obey the Poincare-Hopf relation; 
fN-fB+fR-fc—^ where the N, B, R, and C indices stand 
for the nuclear, bond, ring, and cage critical points [T6j . 



III. STRUCTURAL PROPERTIES 

In order to find the most stable structures of the 
(ZnO)„ (n=2-16) clusters, we calculated and compared 
the optimized energy of several structural isomers of each 
cluster, after applying full geometry relaxation. The ob- 
tained stable isomers are presented in Fig.[T]and their cal- 
culated geometrical properties are listed in Table. HI It is 
seen that the clusters with 9 ZnO units or less prefer two 
dimensional (2D) ring shape geometries while the larger 
clusters stabilize in 3D structures. According to the Ge- 
netic Algorithm the larger clusters are expected to 
be composed of the smaller stable clusters. Our results 
confirm the genetic mechanism as the simple and sym- 
metric rhombus and hexagon structures observed in the 
smaller ZnO clusters are the genetic basis of the stable 
isomers of the (ZnO)9_i6 clusters (Table. H]). In the high 
symmetry (ZnO) 12 and (ZnO)i6 clusters, in spite of their 
different number of atoms, each rhombus is symmetri- 
cally surrounded by four hexagons. In the 11-16 units 
clusters, it is seen that addition of ZnO units to the clus- 
ters mainly increases the number of the hexagon facets 
of the system while the number of the square (rhombus) 
facets is less sensitive to the cluster size. As it was men- 
tioned in the introduction, it is theoretically observed 
that the high symmetry (ZnO) 12 cluster is the building 
block of large ZnO nanostructures [ll| . 

Our results indicate a ring to cage structural cross over 
is between the (ZnO)g and (ZnO) 10 clusters. In fact, the 
(ZnO)9 cluster is structurally distorted from a perfect 
planner geometry (Dgh symmetry) to a zigzag ring with 
Ci symmetry. Hence, the 9 units system is the actual 
onset of the 2D-3D structural cross over in ZnO clus- 
ters. According to the following discussions, the struc- 
tural zigzag distortion of (ZnO)9 is attributed to the be- 
havior of the Zn-O-Zn bond angle (6*0). Because of the 
ionic Zn-0 bonding, the valence electronic shell is mainly 
localized around the oxygen nuclei while the zinc nuclei 
act as the localized positive charge centers in the sys- 
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FIG. 1. The most stable isomers of the (ZnO)n=2-i6 clusters. 
The green (light) and red (dark) balls indicate the Zn and O 
atoms, respectively. The numbers and words written in the 
figure, show n and the point symmetry of the isomers. 



tems. The presence of two lone electron pairs in the va- 
lence shell of oxygen prevents linear configuration of the 
Zn-O-Zn chains (similar to the water molecule) while the 
0-Zn-O chains prefer a linear geometry {9z = 180°). In 
order to estimate the equilibrium value of the angle in 
the ring ZnO clusters, we calculated a free (Zn-0-Zn)^+ 
ligand. The selected charge state of the ligand is based on 
the fact that the nominal atomic ionization in ZnO clus- 
ters is Zn^"'"0^~. Following this method, the equilibrium 
value of 60 is estimated to be about 135°. Inspecting 
the presented data in Table. U clarifies that by increas- 
ing the cluster size, both 9o and 9z angles enhance from 
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TABLE I. Calculated geometrical properties of the most sta- 
ble (ZnO)„ clusters, Sym.: point symmetry, d(A): average 
bond length, 60 and 6z (degree): average Zn-O-Zn and O-Zn- 
O bond angles, frhb and fhex'- number of the genetic rhombus 
and hexagon structural units in the cluster, and /s: number 
of the bond points in the system. 
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small values tovifard their equilibrium values. The 6z an- 
gle increases faster and reaches to about equilibrium lin- 
ear geometry in the (ZnO)5 cluster, while continues to 
increase up to the cluster size of 8. In the perfect planner 
(ZnO)9 structure, is extrapolated to be about 135°, 
close to the estimated equilibrium value. The observed 
zigzag structural distortion in this cluster reduces to 
about 124° (Table, m, hence the real equilibrium value of 
the angle should be about 124°, slightly less than the 
estimated value. In order to verify this equilibrium value, 
we calculated some larger clusters {n = 10, 11, 12) in the 
metastable ring structures and found that all of these 
systems involve the zigzag structural distortion with the 
9o angles of about 124°. The slight difference between 
the equilibrium value of the 9o angle in the 2D ring struc- 
tures and the angle of the free (Zn-0-Zn)^+ ligand may 
be attributed to the closed boundary condition imposed 
in the ring structures. 

The obtained size of the ring symmetry breaking in 
ZnO clusters (n = 9) is considerably larger than that of 
ZnS (n — 5) and ZnSe (n = 4) clusters, while ZnTe clus- 
ters hardly stabilize in the 2D ring shape geometries [l^ . 
These differences are well explained by estimating the 
equilibrium Zn-anion-Zn angles in the 2D ring structures 
of these clusters by using the free ligand calculations. Fol- 
lowing this scheme the Zn-S-Zn, Zn-Se-Zn, and Zn-Te-Zn 
angles were estimated to be about 107°, 104°, and 100°, 
respectively. Since, these equilibrium angles are consid- 
erably smaller than the estimated equilibrium Zn-O-Zn 
angle (135°), the 3D structural distortion in ZnS, ZnSe, 
and ZnTe clusters are expected to occur in smaller sizes 
compared with the ZnO clusters. The difference between 
the equilibrium values of various Zn-anion-Zn angles is 
explained by the Valence Shell Electron Pair Repulsion 
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FIG. 2. Calculated average bond length d, binding energy 
Eb, number of bond points Nb, and single bond energy es 
(= Eb * u/Nb) of the most stable ring (blue circle symbols) 
and cage (red square symbols) isomers of ZnO clusters around 
the ring to cage structural cross over. The vertical dashed 
lines show the position of the cross over. 



(VSEPR) model [20[. In this model, increasing the an- 
ion radius lowers the nuclear attraction on the lone elec- 
tron pairs and consequently increases the lone pairs spa- 
tial distribution, giving rise to a more acute Zn-anion-Zn 
bond angle. 

For more understanding of the ring to cage structural 
transition in ZnO clusters, we have calculated and com- 
pared the structural properties of the metastable cage 
structures of the (ZnO)7_g clusters and the metastable 
ring structures of the (ZnO) 10- 12 clusters (Fig. It 
is observed that in all clusters, the ring structures have 
lower bond lengths and single bond energies and con- 
sequently stronger single bonds, compared with the cage 
structures. The reason is that in the 3D structures atoms 
are more coordinated and hence the valence electrons are 
divided into more bonds. As a result of that, each bond 
receives less electrons and gets weaker. Although, the 
ring structures have stronger individual bonds, the cage 
structures involve more number of bonds. In the small 
sizes, where the number of bonds is very low, the indi- 
vidual bond energy determines the stable structure of the 
system, while increasing the cluster size enhances the role 
of the coordination number and at some point (n = 10) a 
ring to cage structural transition occurs in the system. It 
should be noted that the first 3D cluster (ZnO 10), favors 
an intermediate atomic configuration between the ring 
and cage structures which is composed of two (ZnO)5 
ring structures (Fig.[T]). 

The structural behavior of ZnO clusters are more in- 
vestigated by calculating and studying the topology of 
electron density in these systems. The calculated num- 
ber of the critical points as a function of the cluster size 
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FIG. 3. Right: Configuration of the Nuclear (N), Bond (B), 
Ring (R), and Cage (C) critical points in the (ZnO)8 and 
(ZnO)i2 clusters. Note that the clusters are sketched in dif- 
ferent scales. Left: Calculated number of the critical points 
of ZnO clusters as a function of the cluster size. 

is displayed in Fig. |3] along with the prototype configu- 
ration of the critical points in the (ZnO)8 and (ZnO)i2 
clusters. We observe that the number of the nuclear and 
bond points in the ring clusters (n < 9) are the same, 
because in the ring geometry all atoms have coordina- 
tion of two and hence every bond is shared between two 
atoms. Moreover, as it is expected, the ring clusters have 
only one ring point without any cage critical point. The 
3D clusters {n > 9) have cage like structures (Fig.[T]) and 
hence involve one cage critical point. Based on the num- 
ber of CPs (Fig. [3]) one may distinguish three structural 
regions. The first region belongs to the ring ZnO clus- 
ters {n — 2 — 9) where the number of the bond CPs is 
equal to the number of the nuclei {fs/fN = !)■ The ZnO 
clusters with n = 11 — 14 make the second region where 
all atoms are three coordinated and hence fs/fN = 3/2. 
The (ZnO) 10 cluster, has a transition state between the 
two and three coordination clusters. Transition to the 
third structural region starts at n = 15 where the fs/ In 
ratio exceeds 3/2 toward ratio 2 which corresponds to the 
atomic coordination of four, equivalent to the coordina- 
tion of atoms in bulk ZnO. Therefore it is expected that 
the third region belongs to the bulk like structures or 
some metastable structures (like the sodalite structure) 
predicted for large ZnO nanostructures. 

The obtained binding energies of the most stable 
(ZnO)„ clusters as a function of n are plotted in Fig. H) 
It is recently argued that while the binding energy of 
metallic clusters linearly scales with the inverse of the 
cube root of the cluster size (n~^/^) [2l|, in semiconduc- 
tor clusters, the binding energy decreases much faster 
[2^ . After accurate curve fitting, we realized that the 
binding energy of ZnO clusters obeys the equation: 

EB[eY/fu] = 10.2171-2/3 _ 7.79„-i/3 _ 5 06 
which is a quadratic function of n^^^^. Similar ob- 
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FIG. 4. Calculated binding energy values {Eb) and energy 
second-order differences (A2-E) of the most stable isomers of 
ZnO clusters. The abbreviation fu stands for formula unit. 
The solid line shows the nonlinear curve fitted to the binding 
energies (see text). 

served quadratic behaviors in some other oxide semi- 
conductor nanoclusters including CuO and MgO [2^, 
are attributed to the nontrivial contributions of the sec- 
ond nearest neighbor and higher order interactions in the 
cluster energies. We observe that the binding energy of 
the ring (ZnO)„ clusters (rt = 2 — 9) exhibits faster non- 
linear decay, compared with the cage clusters (Fig. |4]). 
It may be due to the lower coordination of atoms in the 
ring clusters which enhances the impact of the higher or- 
der interactions in the binding energies. Further point, 
the nonlinear curve fitted to the ZnO cluster binding en- 
ergies converges to -5.06 eY/fu for very large clusters 
(n — i> 00), which is higher than the binding energy of 
bulk ZnO (-7.45 eV/fu) Q. This limit may be assigned 
to the sodalite structure proposed as a metastable geom- 
etry for large ZnO nanostructures [ll|. 

In order to compare the relative stability of the (ZnO)„ 
clusters and identify the corresponding magic numbers, 
we followed the conventional procedure of the second- 
order differentiation of the cluster energies. The calcu- 
lated second-order differences, plotted in Fig. 01 show a 
strong peak at n = 12 which indicates the highest stabil- 
ity of the (ZnO) 12 cluster among the studied systems and 
proposes 12 as the lowest magic number of this system. 
Although this finding agrees with the previous theoreti- 
cal studies j23'], experimental measurements state 9, 11, 
and 15 as the lowest magic numbers of ZnO clusters. We 
will demonstrate that, the finite temperature vibrational 
excitations may partly explain this discrepancy. 

IV. ELECTRONIC PROPERTIES 

The calculated HOMO-LUMO (HL) gaps of the most 
stable (ZnO)„ clusters calculated within BLYP are pre- 
sented in Fig. [5l along with the available experimental 
and TDDFT data. We observe that the ring clusters 
(n = 3 — 9) have HL gap of about 2.8 eV while the cage 
clusters {n — 11 — 16) exhibit lower gap of about 2.0 eV. 
Therefore, it is understood that the HL gap is mainly 
sensitive to the cluster geometry, not the cluster size. 
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FIG. 5. Calculated HOMO-LUMO gap of the stable iso- 
mers of ZnO clusters within BLYP before and after apply- 
ing the GW correction, along with the experimental Q and 
TDDFT [I data for comparison. 



and the observed red shift around n = 10 is originated 
from the ring to cage structural cross over in ZnO clus- 
ters. Within the TDDFT scheme, the red shift is found 
to occur in lower sizes (around n = 6). The experimen- 
tal results also indicate a red shift between the (ZnO)5 
and (ZnO) 12 clusters, but due to the lack of any HL gap 
measurement in this range the exact size of the red shift 
and consequently structural cross over is unclear. Al- 
though, the (ZnO)2 cluster has a ring shape structure, 
since its bond lengths and bond angles (Table. HI are far 
from the corresponding equilibrium values, its HL gap is 
considerably lower than the other ring clusters. 

It is observed that the BLYP functional highly under- 
estimates the HL gap values compared with the mea- 
sured data. It is usually attributed to the localized na- 
ture of the conventional exchange-correlation functionals 
which prevents them to correctly describe the excited 
state spectra of materials. The many body perturbation 
based GW approximation is a recent approach proposed 
for improving the excited state properties of molecules 
and semiconductors. In the GW approximation, the 
Dyson expansion of the full green function of the system 
is applied for better treatment of the electronic correla- 
tion. In this scheme, the independent Kohn-Sham quasi- 
particles are allowed to weakly interact via an screened 
coulomb potential and hence the quasi-particle spectra is 
effectively improved toward the direct and inverse pho- 
toemission spectra (PES) [13] ■ For instance the band 
gap of bulk ZnO in the wurtzite structure within GGA 
is about 0.8 eV, significantly less than the experimen- 
tal value of about 3.4 eV, measured by PES [2J|, while 
applying the GW corrections improves the band gap to 
about 4.1 eV f25|. 

We applied this many body based correction to the ob- 
tained BLYP electronic structures and observed (Fig. [5]) 
effective enhancement of the HL gap values without con- 
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FIG. 6. Calculated IR spectra of (ZnO),i clusters in the log- 
arithmic scale. The red starred peaks show the dominant IR 
frequencies. 



siderable modification of their trend. Although, the GW 
corrected gaps are somewhat closer to experiment, they 
are unexpectedly overestimated, while TDDFT seems to 
be more accurate for describing the electronic structure 
of ZnO clusters. In order to understand these differ- 
ences, one should more accurately consider the exper- 
imental techniques used for measuring the HL gap of 
ZnO clusters. The reported experimental data are based 
on the UV absorption spectroscopy which is expected to 
measure the fundamental gap, better described by the 
TDDFT computations, while the GW technique is more 
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appropriate for calculating the real gap of materials mea- 
sured by the photoemission techniques (26)] . Our GW cor- 
rected HL gap values may be comparable to the future 
photoemission measurements on ZnO clusters. The cal- 
culated HL gaps after the GW correction (Fig. [5]) exhibit 
a local maximum in the magic number 12, evidencing the 
high chemical stability of the (ZnO)i2 cluster. 

V. VIBRATIONAL PROPERTIES 

In order to study the vibrational properties of ZnO 
clusters, we applied a displacement of 10~''A to all atomic 
positions of the fully relaxed structures and then calcu- 
lated the forces acting on the displaced atoms and con- 
sequently obtained the Hessian matrices of the systems. 
The vibrational modes were calculated by the diagonal- 
ization of the Hessian matrices and the corresponding IR 
intensities were determined by derivation of the dipole 
moments along these modes. The calculated IR spectra 
of ZnO clusters are shown in Fig. |51 Observation of no 
imaginary frequency indicates the dynamical stability of 
the lowest energy isomers of the (ZnO)„ clusters. We 
observe that the dominant IR frequency (the frequency 
with the highest IR intensity) of the clusters increases 
from n = 2 to n = 9 and then starts to decrease. It might 
be attributed to the observed ring to cage cross over in 
the clusters. The dominant IR frequency of (ZnO)i2 is 
found to be about 546 cm~^, close to the IR frequency 
of diatomic ZnO molecule (537 cm~^) and in agreement 
with recent theoretical results [loj . The experimental vi- 
brational frequencies of bulk ZnO in the ground state 
wurtzite structure are 100, 438 and 584 cm~^ [27| . 

The influence of the thermal vibrational excitations on 
the relative stability of ZnO clusters at elevated temper- 
atures are studied by calculating the vibrational contri- 
butions to the Helmholtz free energy (F) of the systems 
from their obtained vibrational spectra, as follows: 

3N 3N 

F^E+-Y,e, + kBTY, ln(l - e"'^^-) (1) 

i i 

where E is the total energy of the system, i runs over the 
number of vibrational modes, Si is the ith mode energy, 
N is number of atoms in the cluster, ks is the Boltzmann 
constant, T is the kelvin temperature, and /3 equals to 
l/ksT. After obtaining the temperature profile of the 
Helmholtz free energy, we calculated the second order 
differences of this free energy at different temperatures 
and plotted the results in Fig. [71 It is clearly visible that 
while n — 12 is the main low temperature magic num- 
ber of ZnO clusters, increasing temperature enhances the 
relative stability of (ZnO)9 cluster to make it a magic 
system at temperatures above 170 K. As it was men- 
tioned, experimental observations suggest 9,11, and 15 
as the lowest magic numbers of ZnO clusters. Our calcu- 
lated vibrational free energies explain the high stability 
of n = 9 while the magic number 11 is not confirmed 
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FIG. 7. Second order difference of the Helmholtz free energy 
of ZnO clusters as a function of the cluster size at different 
temperatures. The horizontal lines show the reference lines 
at different temperatures. 
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FIG. 8. Calculated heat capacity of ZnO clusters as a function 
of temperature at different sizes (left) and as a function of size 
at two specific temperatures (right). 



by our first-principles calculations. Although, the pos- 
itive value of AF in (ZnO) 15 partly argues stability of 
this system, complete verification of the magic number 
15 needs to the further calculation of larger clusters. 

The obtained Helmholtz free energies are then used to 
calculate the heat capacity of ZnO clusters as a function 
of temperature (Fig. [8]). It is observed that at sufficiently 
high temperatures (T > 170 K), the heat capacity mono- 
tonically increases with the cluster size, while at lower 
temperatures some planner clusters have higher heat ca- 
pacities than some 3D clusters. For more clarification of 
this point, we have extracted the heat capacities of all 
clusters at the two specific temperatures of 50 and 300 K 
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and plotted the results as a function of the cluster size in 
the same figure. It is seen that at T=50 K, the heat ca- 
pacity increases by increasing the size of the ring clusters 
and then drops at the onset of the 2D-3D structural tran- 
sition {n — 9), while at T=300 K a monotonic behavior is 
observed for the heat capacity. We discuss these behav- 
iors by pointing the fact that although the 3D cage clus- 
ters involve more vibrational modes (per atom), the plan- 
ner ring clusters (specially (ZnO)8 and (ZnO)9), because 
of the lower coordination and less bonding of their atoms, 
have more low energy vibrational modes (Fig. 15]). There- 
fore at low temperatures, the ring clusters have more 
active vibrational modes and consequently higher heat 
capacities while by increasing temperature, the number 
of the active vibrational modes of the 3D clusters increase 
faster and hence at temperatures above about 170 K, the 
heat capacities of the cage clusters become larger than 
the ring structures. These observations may also explain 
the high temperature magic number 9 of ZnO clusters. 
While at low temperatures, the cage (ZnO)io cluster has 
a lower energy than the ring (ZnO)g cluster, the thermal 
excitation of high energy vibrational modes reduces the 
relative stability of (ZnO)io and hence makes (ZnO)9 a 
magic cluster at temperatures above 170 K. 

VI. CONCLUSION 

In this paper we applied full-potential numeric atom- 
centered orbital method in the framework of Kohn-Sham 



density functional theory to investigate structural, elec- 
tronic, and vibrational properties of the (ZnO)2-i6 clus- 
ters. It was observed that the most stable isomers of the 
(ZnO) 11-16 clusters are composed of the rhombus struc- 
ture of (ZnO)4 and the hexagon structure of (ZnO)6. We 
found that the clusters with 9 or less ZnO units stabi- 
lize in the 2D ring like structures while larger clusters 
favor 3D cage like geometries. The Zn-O-Zn bond angle 
in the ring (ZnO)9 cluster exceeds its equilibrium value 
and hence induces a 3D zigzag structural distortion in 
the system. The stability of the ring structures in the 
smaller ZnO clusters is attributed to their larger indi- 
vidual bond strength while higher number of bonds in 
the larger clusters enhances the stability of their cage 
structures. Accurate electronic structure computations 
indicate that the ring to cage structural cross over in 
ZnO clusters is accompanied by a red shift in the sys- 
tems. The lowest room temperature magic numbers of 
ZnO clusters were found to be 9 an 12. We argued that 
the high relative stability of the (ZnO)9 cluster at room 
temperature is due to the thermal activation of the high 
energy vibrational modes of the larger clusters. 
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